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1. Introduction 

The determination of properties of functions which are equivariant or invariant under 
the transformations of a compact hnear group (hereafter abbreviated in CLG) G is 
often a basic problem to solve in many physical contexts. G-invariant functions play an 
important role, for instance, in the determination of patterns of spontaneous symmetry 
breaking and structural phase transitions (Landau theory IH 12 El El), in equivariant 
bifurcation theory (see, for instance, ^ and references therein), in crystal field theory 
and in most areas of solid state theory. 

It will not be essentially restrictive, in the following, to assume that G is a 
matrix subgroup of the real group 0(n). In fact, complex linear groups can always 
be transformed into real linear groups through a process of realification and compact 
real linear groups are equivalent to orthogonal linear groups. 

An approach to the study of invariant functions, that fully exploits invariance and 
possible regularity properties, takes advantage of the fact that a G-invariant function 
/, defined on R"-, is a constant along each orbit of G and can, therefore, be considered 
as a function / in the orbit space MP/G of the action of G in M". Geometric invariant 
theory [HIIZI suggests how to get, in principle, / and M"/G. 

The orbit space of a compact group can be realized as a connected semi-algebraic 
subset {i.e., a subset defined by algebraic equalities and inequalities) of a Euclidean 
space M.'^. It turns out to be formed by the union of connected semialgebraic manifolds 
of different dimensions {primary strata). The images of G-orbits with the same orbit- 
type form isotropy-type strata, whose connected components are primary strata. 

The properties of a function defined on ]R"/G, may critically depend on the 
geometry of this space, that has to be known explicitly. A simple way to obtain a 
substantial determination of the algebraic equations and inequalities defining the range 
of an orbit map and its strata has been suggested in [Hlinilin] (see also [HUE])- These 
relations can be expressed in the form of positivity and rank conditions of a matrix P{p), 
whose elements are polynomial functions of p G M''. In this way, the defining relations 
of the strata are obtained in the form of algebraic equations and inequalities. 

In applications or, simply, to get a better understanding of the geometry of the 
stratum (its connectivity properties, its boundary, etc.), one often needs to solve 
explicitly these relations and, sometimes, this is very difficult to do using standard 
algorithms. In this paper we propose a general method to derive the equations defining 
the stratification of orbit spaces of CLG's, in the form of explicit rational parametric 
relations. 

This paper is organized as follows. In section 2 we summarize the results in the 
P-matrix approach to the characterization of the orbit spaces of CLG's jHUH]- In section 
3 we present and prove our results, and in section 4, we illustrate them in a physically 
significant example, related to spontaneous polarization in chiral biaxial liquid crystals. 
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2. An overview of the geometry of linear group actions 

In this section, we shall recall, without proofs, some results concerning invariant theory 
and the geometry of orbit spaces of CLG's (see, for instance, O CH CS] and references 
therein) . 

For our purposes, it will not be restrictive to assume that G is a matrix subgroup of 
the group 0(n), acting linearly in the Euclidean space M" and thus defining the G-space 



We shall denote hy g ■ x the action of (7 G G on a; G M" and by the isotropy 
subgroup of G at X G M". The isotropy subgroups of G at points of the orbit G ■ x 
form the conjugacy class (G^) of G^ in G, which identifies the orbit-type of (the points 
of) G • X. The set of points x G (of G-orbits) with the same orbit-type forms an 
isotropy-type stratum o/(G, M""). 

The orbit space of the action of G in is defined as the quotient space M'^/G, 
endowed with the quotient topology and differentiable structure. The image of a stratum 
of M", through the canonical projection vr : M" — > M"/G, defines an isotropy-type 
stratum of M"/ G; the connected components of an isotropy-type stratum of M"/ G are iso- 
dimensional manifolds {primary strata of M^/G), but the orbit space is not a manifold. 

Almost all the points of M"/G belong to a unique stratum, the principal stratum, 
which is a connected open dense subset of M"/ G. The boundary of the principal stratum 
is the union of disjoint singular strata. Every stratum S of M"/ G is open in its topological 



The following partial ordering can be introduced in the set of all the orbit-types: 
(H) < (K), if H is conjugate with a subgroup of K. If (Hi) is the orbit-type of the 
stratum Si of R"/G , z = 1, 2, then Hi < H2 iff S2 is in the boundary of 5*1; therefore, 
more peripheral strata of ]R"/G have greater orbit-type. The number of distinct orbit- 
types of G is finite and there is a unique minimum orbit-type, the principal orbit-type, 
corresponding to the principal stratum. 

The ring R[a;]*^ of real G-invariant polynomial functions of x is finitely generated 
El CHICS]- There exists, therefore, a finite minimal collection of homogeneous G- 
invariant polynomials p{x) = {pi{x) , p2{x) , . . . , Pg{x)) (minimal integrity basis for the 
ring of G-invariant polynomials, henceforth abbreviated in MIB) such that any element 
F G can be expressed as a polynomial function F of p{x): 



The number q of elements of a minimal integrity basis and their homogeneity degrees 
(ij's are only determined by the group G. 

The group G is said to be coregular if the elements of its MIB's are algebraically 
independent. The elements of a MIB of a non coregular group satisfy a set of algebraic 
identities in R": Fa{p{x)) = 0, A = 1, . . . and the associated set of equations 



(G,R'^). 



closure S. 



F{p{x)) = F{x), Vx G W 



(1) 



Fa{p) = 0, 



A 



1,... 



(2) 
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defines an algebraic variety in R'', which is called the variety Z of the relations (among 
the elements of the MIB). If G is coregular, Z = W. 

Since G is a compact group, the orbits of G are separated by the elements of a 
MIB of G (for an elementary proof see, for instance, j^), i.e., at least one element of a 
MIB of G takes on different values at two distinct orbits. Thus, the elements of a MIB 
of G yield a parametrization of the points of MJ^/G, that turns out to be also smooth, 
since the orbit map p : M" — > M'^, which maps all the points of M" lying on an orbit 
of G onto a single point of M'^, induces a diffeomorphism of MT'/G onto a semialgebraic 
g-dimensional connected closed subset of M"^. Like all semialgebraic varieties PU^, p(M") 
presents a natural stratification in connected semialgebraic sub-varieties, called primary 
strata, which turn out to be the connected components of the isotropy-type strata^. 

A characterization of the image p(]R"') of the orbit space of G as a semi-algebraic 
variety and of its can be easily obtained through a matrix P{p), defined 

only in terms of the G-invariant Euclidean scalar products between the gradients of the 
elements of the MIB {p{x)}: 

n 

Pabix) = ^ diPaix) diPb{x) = Pabipix)), a,b = l,...,q, (3) 

1=1 

where in the last member, use has been made of Hilbert's theorem, in order to express 
Pafc(x) as a polynomial function of pi{x), . . . ,Pq{x). 

The following theorem [HI El CHI (see also clarifies the meaning and points out 
the role of the matrix P{p): 

Theorem 2.1 Let G be a compact subgroup of the real group 0(n), p the orbit map 
defined by the MIB {pi{x) , p2{x) , . . . ,pq{x)) and P{p) the matrix defined in 
([^. Then the range p{W^) of p is the unique semialgebraic connected subset of the 
variety Z <ZW^ of the relations among the elements of the MIB where P{p) is positive 
semi-definite. The k- dimensional primary strata ofp{W^) are the connected components 
of the set {p E Z \ P{p) > 0, rank(P(p)) = k}; they are the images of the connected 
components of the k-dimensional isotropy-type strata ofW^/G. In particular, the set of 
interior points o/p(M"), where P{p) has the maximum rank, is the image of the principal 
stratum. 

In the following, we shall identify orbit spaces and their strata with their images 
through orbit maps. 

If G is coregular, the Pi's, with range in the semialgebraic set A defined by the 
inequalities P{p) > (and no equation!) provide a one-to-one parametrization of the 
principal stratum of M"/G. If the stratum is singular, the pi's are not independent 
parameters, being bounded also by the equations defining the stratum. In applications 

I A simple example of a compact connected linear semialgebraic variety of is yielded by a 
polyhedron: its interior points form a unique three-dimensional primary stratum (principal stratum), 
while two-, one- and zero-dimensional primary strata are formed, respectively, by the interior points 
of each face, by the interior points of each edge, by each vertex. 
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or, simply, to get a better understanding of the geometry of the stratum (its connectivity 
properties, its boundary, etc. ...), one often needs to solve explicitly these equations, 
which arc obtained in implicit form from the theorem recalled above. Often, this is very 
difficult to do using standard algorithms. In the following section, we propose a general 
method to derive the equations defining the singular strata of the orbit space of any 
CLG in the form of explicit rational parametric relations. 

The mathematical apparatus just recalled offers effective tools, for instance, in the 
determination of possible patterns of spontaneous symmetry breaking, when the ground 
state of the system is determined by the absolute minimum of an invariant potential. 
Let us recall the basic ideas. In this context, the vector x e is an order parameter 
and G is the symmetry group of the potential ^{a;x) (free energy, or Higgs potential, 
for instance), expressed in terms, also, of parameters a. The points XQ{a), where the 
function (f)a{x) = ^{a;x) takes on its absolute minimum, determine the stable phase of 
the system, whose residual symmetry is defined by the isotropy subgroup of G at Xq. 
Owing to the G-invariance of the potential, each of its stationary points is degenerate 
along the G-orbit through it. Since the isotropy subgroups of G at points of the same 
orbit are conjugate in G, only the orbit-type of xo{a) is physically relevant. Structural 
phase transitions take place when, by varying the values of the a's, the point Xo{a) is 
shifted to a different stratum. 

If $(a; x) is a sufficiently general function of the a's, by varying these parameters, 
the point .To(a) can be shifted to any stratum of M^/G. So, the strata are in a one-to-one 
correspondence with the symmetry phases allowed by the G-invariance of the potential. 
On the contrary, extra restrictions on the form of the potential function, not coming 
from G-symmetry requirements (e.g., the assumption that the potential is a polynomial 
of low degree), can limit the number of allowed structural phases for the system in its 
ground state. 

3. Parametrization of strata in orbit spaces of compact linear groups 

Let G be a compact linear group, acting orthogonally in and {p} be a related MIB. 
We shall prove that the possibility of parametrizing the principal stratum of the orbit 
space of a coregular group in terms of parameters related to the elements of a MIB can be 
extended to singular strata. The result stems from the proof of the following statement: 
the set of interior points of the topological closure of an isotropy-type stratum, with orbit- 
type (H), is diffeomorphic to the principal stratum of the orbit space of a group space of 
the stabilizer of H inG. If the stabilizer turns out to be coregular. the constructive proof 
of this proposition will provide a one-to-one rational parametrization of the stratum. 
Since, in view of possible applications, we are mainly interested in the possibility of 
getting a one-to-one parametrization of a stratum, in the presentation of our results, 
stress will be laid on this aspect of the matter. 

Let if be a proper isotropy subgroup of the G-space (G, W). We shall denote by S 
the stratum of (G, M") with orbit-type {H) and by S its image p{S) in the orbit space 
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MJ^/G. Let us, also, define 

V = {x eW \G^^ H}, V = {xeW'\G^ = H}. (4) 

The non trivial sets V and V have the following properties, which are more or less 
immediate consequences of their definitions: 

(i) V = {x E MJ^ \ h ■ X = X, h E H} is a linear subspace of M", let us call u its 
dimensions; 

(ii) V is the topological closure of V: V = V; 

(iii) V = S nV and every G-orbit lying in S has at least a point in V so that S = G - V, 
S = G ■ V (where the bar denotes topological closure) and, consequently, 

p{V) = S, p{V) = S (5) 
and 

rank(P(a;)|^g^) = dim(5), Ta.iak{P{x)\^^y^^) < dim{S). (6) 

A G-orbit of 5* may intersect V in one or more points. In the first case, every G-orbit 
of S intersects V in a point and, owing to item (iii), the coordinates of the points of V 
provide a one-to-one rational (in effect, linear) parametrizazion of the orbits of G lying in 
S, obtained by solving the system of linear equations h-x = x\fh E H . The allowed range 
of X, X G V, is determined by the inequalities assuring that rank(P(x)) =dim(5'). The 
parameters x, with range V, would provide, in this case, a one-to-one parametrization 
of the stratum, thus solving our problem. In general, however, the intersection of a 
G-orbit, of orbit-type (if), with the stratum S does not reduce to a single point. So, a 
sounder analysis is required, that will be the object of the rest of this section§. 

Two distinct points, x and g ■ x, g G G, of the same G-orbit, lie in V iff 
h ■ g ■ X = g ■ X "^h E H, that is, iff g is in the stabilizer Stab (if, G) of H in G. 
The intersection of a G-orbit of S with V is, therefore, the Stab (if, G)-orbit through x 
and Stab(ii, G) is the largest subgroup of G leaving V invariant. 

In the group space (Stab(if, G),V), the isotropy subgroup at a point of general 
position is H. Therefore, the principal stratum S satisfies the following relation: 

vdj:^v. (7) 

It has to be noted that S could contain V in a strict sense, since, at points of V\V, the 
conjugacy class in G of the isotropy subgroup of Stab(ii, G) could be smaller than the 
conjugacy class in G of the isotropy subgroup of G. 

Now, H is an invariant subgroup of Stab (if, G) and a subgroup of all the isotropy 
subgroups of (Stab(if, G), V). Thus, the action of Stab(if, G) in V defines a linear 

§ Despite its being, generally, not one-to-one, the parametrization of the orbits of a stratum by means 
of the points of V may be useful. It has been used, for instance, by Kim 21 to parametrize the strata 
of a set of low dimensional orbit spaces of compact linear Lie groups. 
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group K (and the group space {K, V)), isomorphic to the quotient group Stab(if, G)/H, 
through the relation 

(sH) ■v = s-v, se Stah{H, G), sH e Stab(/7, G)/H. (8) 

So, we can conclude that the intersection of a G-orhit of S with V is an orbit of {K, V). 
Obviously, the orbit spaces V/K and V/Stab{H, G) are isomorphic and can be identified. 

We can rephrase the result just obtained, by claiming that the points of S are 
in a one-to-one correspondence with the points of a, possibly proper, subset V of the 
principal stratum S of the orbit space V/ K and the following relations hold true: 

S = p{V) D D p{V) = S. (9) 

Since, as stressed, we know how to parametrize a principal stratum, our problem is 
reduced to the determination of V. We shall show that this can be easily done, making 
use of dni). 

To make simpler the solution of the problem, let us assume that an orthonormal 
basis has been introduced in such that the first u elements of the basis yield a basis 
for the vector space V. Then, if we denote by V± the orthogonal complement of V in 
M", the subspace V± is invariant under [K and) H, owing to the orthogonality of the 
transformations of G. Since V contains all the if-invariant vectors of M", there is no 
non trivial if-invariant vector in V±. 

To attain our goal, we shall start from a convenient parametrization of the principal 
stratum S of V/K in terms of / real parameters A, related to a MIB (Ai(f ),..., A/ (f)) 
of the ring of polynomial /^-invariant functions of v & V. This parametrization will be 
global and one-to-one, if K turns out to be coregular. In this case, the range of A has 
to be restricted to the positivity region of the P-matrix A(A), associated with the MIB 
{A}: 

V 

Aa/3(A(t;)) = Ao/3(t^) = ^ <9iA«(a;i, ...,Xy) diXp{xi, ...,x^). (10) 

1=1 

As recalled in the previous section, in fact, the orbit space V/ K and its principal stratum 
E can be identified, respectively, with the semialgebraic sets XiV) and A(S): 

S = A(S) = {A G M' I A(A) > 0}, V/K = S = X{V) = {XeR^\ A(A) > 0}, (11) 

the second set being the closure of the first one. Moreover, the definition of V and 
equation imply 

V = A(V) C A(S) = S. (12) 

If K is not coregular, only a local one-to-one parametrization can be obtained for 
S, by eliminating redundant elements in the set of parameters (Ai, . . . , A;), through the 
solution of the algebraic relation(s) among the elements of the MIB {A}, and imposing 
convenient semi-positivity and rank conditions on the matrix A(A). 
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A one-to-one local or global parametrization of S yields, obviously, also a local or 
global one-to-one parametrization of V, provided that additional restrictions are imposed 
on the range of A, whenever V is a proper subset of E. In this case, the correct bounds 
can be obtained in the following way (we shall only consider the case of a coregular K, 
the extension of the results to non coregular K^s is straightforward, but, as just stressed, 
may lead to a loss of globality). 

When X spans V, the elements of the MIB {pi{x), . . . ,Pq{x)) of G define a set of 
i^-invariant polynomial functions of f = {xi, . . . ,Xu) G V. Therefore, by the Hilbert 
theorem recalled in the previous section, they can be expressed as polynomial functions 
of A, that is, 

p\^ = (f)oX. (13) 

Let us remark that, in our assumptions, there are no relations among the elements of the 
MIB {A} and, consequently, possible relations Fa{p) = among the pj's are identically 
satisfied for p = 0(A). 

From p3|) and (|7j), one immediately obtains 

S = p{V) = (j){V/K) D 0(S) D 0(V) = p{V) = S (14) 

and, since S, being a principal stratum, is the set of interior points of V/K and is 
connected, 0(S) will coincide with the set of interior points of the closure S" of and 
will be connected. This set does not coincide with S* if S* contains, in its interior, points 
representing bordering strata of S. This certainly happens if the set S is not connected 
and, presumebly, also if it is not multiply connected||. 

The identification of points A G S, if any, whose image 0(A) ^ S, can be obtained 
in the following way. 

Let X = V (B v± yield the decomposition of x G in its vector components v & V 
and v± G V±: v = (xi, . . . ,x^), v±_ = {x^+i, . . . ,x„). Then, a G-invariant polynomial 
f{x), can be thought of as a polynomial in v and v± and it is easy to realize that it 
cannot contain linear terms in v±, being v invariant under H. As a consequence, 

dif{x) = 0, for X G V and i = u + 1, . . . ,n. (15) 

Now, starting from the very definition of P{x) (see for every x = f G V" we obtain, 
using (fTHjl and the identity p{v) = 0(A(f)), 

Pab{p{v)) = Yli=i diPa{v) diPb{v) 
J(A)A(A) J(A)T 

A=A(ti) 

II To our knowledge, no general proof exists even of simple connectivity of the principal stratum of a 
coregular compact group; however, principal strata of coregular low dimensional (D < 4) orbit spaces 
can be checked to be multiply connected |22l 1281 . 
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where, the superscript T denotes transposition and J(A) is the Jacobian matrix of the 
transformation p = 0(A): 

JaaW = da(t>aW, a = l,...,q, a = l,...,l. (17) 
So, we can conclude that, for all \ &V/K, 

P(0(A)) = J(A)A(A)J(A)T. (18) 

Equation (fTHj) leads to an easy calculation of the points A G V, that is, of the points 
A G S whose image p = 0(A) is in S. In fact, from (jHl), these points are characterized 
by the conditions A G S and rank(P(0(A))) = /. Since, for A G S, the matrix A(A) is 
positive definite with rank /, we can conclude that the range of A assuring the location of 
the point p = 0(A) in S is the semialgebraic subset A of M', determined by the following 
inequalities: 

A = {A G I A(A) > and rank(J(A)) = /}. (19) 

For A G A, the relation p = 0(A) yields a global rational parametrization for 5*. 

It will be worthwhile to stress that the parametrization we have suggested is, in 
some way, canonical: the unique arbitrariness in the choice of the parameters is related 
to the choice of the MIB's. 

An important byproduct of the result just proved is a simple test of the connection 
of 5*, that could be difficult or impossible to read directly from the equations of the 
stratum in implicit form. In fact, the condition that the boundary of S coincides with 
the boundary of S is equivalent to the condition rank(J(A)) = / for all A G S. 

If K is not coregular, there are polynomial relations -Fa (A) = among the elements 
of the MIB (Ai,...,Ai) and / >dim(S'). As already stressed, to obtain a one-to-one 
parametrization of the points of S by means of the A^'s, one has to eliminate the 
redundant parameters by solving the equations -Fa (A) = 0. This may be feasible, 
but only locally and the resulting parametrization will not be global and possibly not 
rational. 



4. An Example 



In this section we shall show how the parametrization technique works in a simple 
example. The notations will be the same defined in the previous sections. 

We shall consider the orthogonal linear group G defined by the action of the real 
group 0(3) in the real eight-dimensional space spanned by the independent components 
Xi, . . . , X5 of a symmetric and traceless tensor Q and the three components Xg, x-r, Xg of 
a polar vector P. To be specific, if O is a generic 3x3 real orthogonal matrix, the 
transformation rules of the Xi are obtained from the following relations: 



/-A 

X3 



Xi 



X3 



X4 



\ 



Xi 

7! 



- X2 



X5 



(20) 



X4 



X5 



+ X2 
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Q' 



p' 

a 



Yli'y5=l ^"7 QlS Ops 
Yl^f3=l Oap P(3 , 



a,(3 = 1,2,3. 



(21) 



It is worthwhile to remark that the reflection / = diag(— 1, —1,-1) G 0(3) reverses only 
the signs of the last three coordinates (xe, x-j, xg) and, due to the symmetry of the tensor 
Q, each G-orbit contains points where Q takes on a diagonal form {x^ = X4 = X5 = 0). 
This remark makes much easier the determination of "typical points" in strata, from 
which the orbit type of the stratum can be identified I21j. 

The real orthogonal linear group G, just defined, is coregular (see j22j) and its 
isotropy subgroup at a generic point of is trivial. As a consequence, the principal 
orbits are three-dimensional manifolds, the principal stratum has dimensions five and 
there are five independent elements in a MIB (see, for instance PHj), which can be 
chosen in the following way: 



Pi = Tr g2 + p . p 



8 




P2 = P P 



8 




P3 = QV2 Tr Q 



,3 



2 VSxi^ + 6 y/Sxi X2 — 3 \Pixx — 9 0:2 x^^ 

3 xi + 9 X2 xa^ + 18 X3 X4 X5 + 6 Xi x^ , 



(22) 



Pi = 3V2^Po,QapPp 



= — 2 -\/3 xi Xq' + 6 X3 Xe X7 + \/3 Xi Xy^ — 3 X2 Xj^ 
+ 6 X4 xe X8 + 6 X5 X7 xs + xi Xg^ + 3 X2 Xg^ , 



P5 = 6j2PaQ\pPp 



= 4 Xi^ Xg^ + 3 X3^ Xq^ + 3 X4^ Xg^ — 2 Xi X3 Xg Xy 

— 6 X2 X3 Xg X7 + 6 X4 X5 Xg X7 + Xi^ Xy^ — 2 -\/3 Xi X2 X7^ 
+ 3 X2^ Xy^ + 3 X3^ X7^ + 3 Xs^ Xy^ — 2 -\/3 Xi X4 Xg Xg 

+ 6 X2 X4 Xg Xg + 6 X3 X5 Xg Xg + 6 X3 X4 X7 Xg + 4 Xi X5 X7 Xg 
+ xi^ xg^ + 2 -\/3 Xi X2 Xg^ + 3 X2^ Xg^ + 3 X4^ Xg^ + 3 Xs^ Xg^ . 
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The corresponding P-matrix elements can be easily calculated from their definition Q: 

Pla = '^■daPa, l<a<5 

P22 = 4p2 

P23 = 

P24 = 4p4 

P25 = 4p5 

P33 = 108(pi-p2)' (23) 

P34 = 18{-2pip2 + 2p2^+p,) 

P35 = 12 (P2P3 -P2P4) 

P44 = 12 fe'+Ps) 

P45 = 4 (P2P3 + 3piP4 -P2P4) 

= I (P3P4 +P4^ + QpiPs) , 

where the ordered set (rfi, (i2, c?3, ^4, ^5) = (2,2,3,3,4) specifies the degrees of the 
polynomials of the MIB. 

The determinant of the matrix P{p) factorizes and only one of the two real 
irreducible factors, that we shall call A{p), turns out to be active 24j, that is, to be 
relevant in the determination of the boundary of the orbit space M.^/G: 

A{p) = 3p2^P3'^ + I8pip2'^ PsPi- 18 P2^PS Pi + 27 Pi'^P2P4^ 

-54piP2^P4^ + 27p2^P4^ +P3P4^ - 9P2P3P4P5 (24) 
-9piP42p5 + 9p2P4^P5 - 27piP2P5^ + 27^2^ Ps^ + 9^5^ . 

The relations defining strata of dimension < 4 in the orbit space of (0(3),]R®) are 
summarized in tablelH The isotropy subgroup lattice, with the possible phase transitions 
between bordering strata is shown in figure H 

4.I. Parametrization of the strata 

The relations assuring that P{p) > and has rank 4, define a unique four-dimensional 
stratum S^'^^ in the orbit space. Using well known matrix theory results, these conditions 
could be explicitly written, for instance, in the form A{p) = and Mi{p) > 0, 
i = 1,...,4, where Mj is the sum of the principal minors of order i of the matrix 
P{p): a cumbersome set of conditions that it is not worthwhile to write down explicitly. 

Even in this simple example one immediately realizes that the choice of a minimal 
set of explicit algebraic relations providing a cylindrical decomposition [23] for the 
semi-algebraic subset S^^^ of would be a really hard task (for the more peripheral 
strata, instead, the problem is much easier to solve). An immediate application of the 
results proved in the previous section, on the contrary, leads to a simple rational global 
parametrization of each stratum, as shown below. 

4- 1.1. Stratum S^'^\ A "typical point" of the stratum is Xt = (1, 1, 0, 0, 0, 0, 1, 1). The 
corresponding isotropy subgroup if of G is the Z2 group generated by the refiection 
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representing the element diag(— 1, 1, 1) G 0(3) in M^: 

diag(l,l,-l,-l,l,-l,l,l)GG. 

The vector space V formed by the if-invariant vectors of M*^ turns out to be five- 
dimensionah 

V = {x eR^ \ X3 = Xi = Xq = 0} . 

The elements of 0(3), corresponding to elements of Stab [H, G), are block- diagonal 
matrices of the form diag(±l, O), with O G 0(2). Therefore, Stab {H, G) is isomorphic 
to a group Z2xO(2) and the quotient group K = Stab {H,G)/H is the representation 
in V (induced by the representation G of 0(3)) of the Oi(2) subgroup of 0(3), formed 
by the rotations around the first axis. 

Using coordinates v = (xi, X2, X5, xy, xs) for a vector v E V, the elements of K 
turn out to be block diagonal matrices leaving invariant the following subspaces V^, 
1<A<3: 

W = {(xi, 0,0, 0,0) G V\xi G M}, 

= {(0, X2, X5, 0, 0) G V\X2 , X5 G M} , 

= {(0, 0, 0, X7, xg) G V\x7 , xg G M} . (25) 

A proper rotation r(0) G Oi(2) of an angle and the reflection diag(-l,l)G Oi(2) 
are represented in V by the 5x5 matrices 1 © r(— 20) © r(0) and diag(l, —1, 1, —1, 1). 

Noting that the complex variables zi = X2 + ixs and Z2 = xr + ixs transform into 
exp(— 2i0) zi and, respectively, exp(i0) Z2 under a rotation and into z^ and, respectively, 
—Z2 under a reflection, it is easy to realize that a possible MIB for {K, V) is the following: 

Ai = xi , 

\ _ I |2 _ 2 I 2 
^2 — 1^1 1 — + X5 , 

\ I |2 2,2 

A3 = F2I = + Xg , 

A4 = 2Re(2;iZ2^) = 2 (x2X7^ — 2x5X7X8 — X2X8^) . (26) 

It is, now, easy to express the p's in term of the A's, p = 0(A), and to check that the 
following expressions, obtained in this way, identically satisfy the equation A{p{X)) = 
(see fl24|l ) entering the definition of the stratum S^^^: 

0i(A) = Ai' + A2 + A3, 

02(A) = A3 , 

03(A) = -2v^Ai (A12-3A2) , 
04(A) = ^ (2A1A3 - V3\,) , 

05(A) = Xi^Xs + 3A2A3 - V3A1A4. (27) 

As explained in the previous section, since the group K is coregular, the range A for 
A is the region where the P-matrix A(A) associated to the MIB {A} is positive definite 
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and the rank of the Jacobian matrix J(A) of the transformation 0(A) is maximum (=4). 
The matrices A(A) and J(A) are easily calculated to be 



A(A) 



and 



J(A) 



/ 1 





\ 











4A2 


2A4 











4A3 


4A4 








Vo 


2A4 4A4 4 (A3 


' + 4A2A3) ) 










2Ai 


1 


1 





\ 










1 









6^3 (Ai'- A2) 


e^sAi 












V^A3 







-3/2 




V 


2A1A3 — \/3A4 


3A3 Ai^ 


+ 3A2 


-v^Ai 


/ 



(28) 



(29) 



The conditions assuring the positivity of A(A) > can be written in the form 
A = {(Ai, A2, A3, A4) G I A2 > and A3 > and Xl < 4X2X3^} 



(30) 



and the rank of J(A) turns out to equal 4 for all A G A. So we can conclude that the 
stratum is connected, a piece of information that would be difficult to derive directly 
from the relations defining the stratum in implicit form {i.e. A{p) = and Mi{p) > 0). 

The parametrization obtained for the sub-peripheral stratum S^^^ turns out to be 
useful, also, because bordering values of A immediately determine bordering values 
of p, corresponding to more peripheral strata of M.^/G. The stratification of V/K is 
summarized in table |21 We would like to remark, however, that there is not a one-to- 
one correspondence between singular strata of {K, V) corresponding to bordering values 
of A and singular strata of M"/ G, corresponding to the associated values of p{X). In fact, 
as was already noted, since the action of the stabilizer on the subspace is trivial, 
the invariant Ai has degree 1, so it does not participate in the conditions defining the 
stratification of {K,V). That is why the parametrization procedure has to be applied 
stratum per stratum. 

4.1.2. Stratum S^^\ A typical point in this stratum is Xt = (1,1,0,0,0,1,0,0) and 
the isotropy subgroup at is the Z2 x Z2 group represented by the 0(3) elements 
diag(l, itl, ±1). The vector space left invariant by H reduces to the three-dimensional 
space V: 

= {x G I X3 = X4 = X5 = X7 = Xg = 0} . (31) 

The linear group acting in V , is a Z2 x Z2 group, generated by the matrices 
diag(l, -1, 1) and diag(l, 1, -1). A MIB for Stab {K, V) is 

(Ai, A2, A3) = (xi, X2^ X6^) . (32) 
The associated matrix A (A) turns out to be 

A(A) = diag(l,4A2,4A3). (33) 
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Parametric equations, p = 0(A), for the stratum can be written in the following form: 

0i(A) = Ai' + Aa + As, 

02(A) = A3 , 

03(A) = -2v^Ai (Ai'-3A2) , 
04(A) = -2V3A1A3, 

05(A) =4Ai2A3, (34) 

for A in the range A = {(Ai, A2, A3) G | A3 > and A2 > 0}. 

The rank of the jacobian matrix J(A) turns out to equal 3 for all A G A, so that 
the stratum is connected. 

4.1.3. Stratum S^^\ A typical point is Xt = (1,0,0,0,0,1,0,0) and the isotropy 
subgroup if of G at Xt is the subgroup formed by the rotations (proper and improper) 
around the first axis, represented by block-diagonal matrices of the form diag(l,0), 
with O G 0(2). The vector space V left invariant by H is two-dimensional: 

= {x G I X2 = X3 = X4 = x^ = X7 = Xs = 0} . 

K is the Z2 group generated by the matrix diag(l,— 1). A MIB for {K.,V) is the 
following: 

(Ai, A2) = (xi,X6^) . (35) 
The associated matrix A(A) turns out to be 

A(A) = diag(l,4A2). (36) 
Parametric equations, p = 0(A), for the stratum can be written in the following form: 



2 : 



0i(A) =Ai^ + A 

02(A) = A2 , 

03(A) = -2v^Al^ 
04(A) = -2V3A1A2, 

05(A) =4Ai%, (37) 

for A in the range A = {(Ai, A2) G | A2 > 0}. The rank of J(A) is equal to 2 on the 
whole A. 

4.1.4. Stratum S]^\ A typical point is Xt = (1,1,0,0,0,0,0,0) and the isotropy 
subgroup if of G at Xt is the subgroup Z2 x Z2 x Z2 formed by the reflections of 
the axes in M^. The vector space left invariant by H reduces to 

V = {x G M.^ \ X3 = Xi = X5 = xq = Xt = xs = 0} . 

The group i^' is the finite group of order 6 generated by reflections with root system of 
type A2 [23^ with generators 

/ -1/2 V3/2 A / -1/2 -VS/2 \ 



[ V3/2 1/2 j ' V3/2 -1/2 ) 
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A MIB for Stab {H, G) is 

(Ai, A2) = (xi^ + X2^ , - 3X2^)) . 

The associated matrix A (A) turns out to be 



(38) 



A(A) 



4Ai 6A2 
6A2 QAi^ 



) 



(39) 



Parametric equations, p = </'(A), for the stratum can be written in the following form: 



for A in the range A = {(Ai, A2) G | Ai > and Ai^ - Aa^ > O}. The rank of J(A) is 
equal to 2 for all A G A. 

4.1.5. Stratum S^^\ A typical point is Xt = (1,0,0,0,0,0,0,0) and the isotropy 
subgroup if of G at Xt consists of the 0(3) block-diagonal matrices of the form 
diag(±l,0), with O G 0(2). The vector space left invariant by H reduces to 

V = [x & \ X2 = X3 = X4 = X5 = Xq = Xt = Xs = 0^ , 

the group K is trivial and a MIB for Stab {K,V) is simply Ai = Xi. The associated 
matrix A(A) reduces to 1. 

Parametric equations, p = 0(A), for the stratum can be written in the following 
form: 



for A in the range A = M. One can verify that the rank of J(A) diminishes in the allowed 
range A for Ai = 0. 

4.2. A physically interesting non coregular variant of the example 

The example presented in the first part of this section has been suggested by a paper of 
Longa and Trebin [2H] , devoted to the construction of a phenomenological theory of polar 
structures in chiral biaxial hquid crystals, through the exploitation of the properties, 
under transformations of a symmetry group SO (3), of a symmetric and traceless tensor 
order parameter field Q and of a polar vector field P. 

In fact, the situation examined by Longa and Trebin can be recovered from the 
example worked about, considering the action in of the subgroup S0(3) of 0(3). 
In this case, we obtain a non-coregular linear group G' . To get a MIB for G", a new 
G'-invariant polynomial p&{x), of degree 6, has to be added to the MIB (j22I): 



(0i(A), 02(A), 03(A), 04(A), 05(A)) = (Ai , , -2 v^A2 , , o) , 



(40) 



(0i(A), 02(A), 03(A), 04(A), 05(A)) = (Ai' , , -2 Al^ , 



') 



(41) 



P6 = 2y2(PAgP)-(Q2p) . 



(42) 
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The P-matrix associated to the enlarged MIB can be constructed from the one shown 
in by adding one more row and column: 

^16 = 12^6, 
P26 = 6p6, 
^36 = 0, 
^46 = 0, 
P56 = 12piP6 , 

-Pee = - 7: [-Pi'^Pi^ +P4 {P3+P4) Pb - ^PiPb^ + 2pip2 {-P3P4 -P4^ + QP1P5) + 

(ps' - 2 p3 P4 - 3 P4' + 12 pi ps) ] • (43) 

The added element ^6(2;) is a "numerator invariant" for the Molien function [21] and 
is, therefore, algebraically dependent on the other elements of the MIB. The relevant 
relation can be easily obtained from one of the irreducible polynomial factors of the 
enlarged P matrix: 

243 pe^ + A{p) = 0, (44) 

where A{p) is the same (see fl24j) ) as in the coregular case. 

Being G' a subgroup of G, the lattice of isotropy types of (S0(3),M^) is easily found 
from the following relation, holding true for all x G M^: 

G', = G,f]G'. (45) 

An explicit calculation shows that, the number of strata of (80(3),^*^) is lower 
than the number of strata of (0(3),M^). In fact, the four-dimensional stratum S^^^ of 
the second group space is part of the principal stratum of (S0(3),M^). The result is 
summarized in table El where the notations are the same used for the generators in the 
0(3) case. 

As explained at the end of section 2, if SO (3) is assumed to be the (largest) 
symmetry group of the free energy, the strata of M®/G" are in a one-to-one 
correspondence with the structural phases of the system, i.e., with the phases of 
the system that are identified only on the basis of their SO (3) symmetry. An easy 
comparison with the results summarized in table I (p. 3460) of jSHj, shows that distinct 
polar states (Pq}^, Fsg and Fu2) defined by Longa and Trebin, lie on a same stratum 
(the principal stratum) of (S0(3),M^), meaning that their symmetry (orbit type) is the 
same: they form a unique structural phase. As stated above, the number of structural 
phases is increased if the symmetry group of the free energy is enlarged to 0(3) (see our 
table E}, but also in this case the polar states Fsg and Fu2 lie on the same stratum S^'^\ 

In their paper, Longa and Trebin make a big effort to determine the range of the 
orbit map and produce a classification of the polar states that, however, appears to 
be not based on their symmetry properties with respect to the symmetry group of the 
free energy. The effectiveness and mathematical rigor of the orbit space approach to the 
study of invariant functions clearly emerges from a comparison with the easy calculation. 
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that led us, essentially, to more rigorous results, and with the additional important 
advantage of keeping, at every step, a clear correspondence between geometrical 
structures (strata) and symmetry of the corresponding physical configurations. In fact, 
the parametrization technique becomes very important for orbit spaces of group actions 
where more than three invariants are involved; here the geometrical intuition based on 
drawing the shape of the orbit space loses most of its efficiency. 
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Table 1. Relations defining strata of dimension < 4 in the orbit space of (0(3),M^). 



Stratum Equalities Inequalities 



12pip2' - 12p2^ > 
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Table 2. Isotropy type stratification of the orbit space V/K, whose principal stratum is difFeomorphic 
to the subprincipal stratum S'*^''-' of the orbit space R^/G. 



Stratum 


Equalities 


Inequalities 


E(2) 


A3 = A4 = 


A2 > 


E(3) 


4A2A3^ = 


A3 > and A2 > 






A2 > and A3 > and 4A2A3^ - A4^ > 



Table 3. Strata and isotropy types of the space (SO (3), M^). For each orbit type, a typical point and 
the corresponding isotropy subgroup (or its generators) are specified, with reference to the form of the 
matrices O appearing in (|21|) : 0± denote elements of 0(2) with determinant ±1, respectively. 



Stratum 


typical point 


isotropy subgroup or set of j 


generators 


^(5) 


(1,1,0,0,0,0,1,1,1) 


{e} 






(1,1,0,0,0,0,1,0,0) 


{diag(l,-l,-l)} 




o{2) 


(1,0,0,0,0,0,1,0,0) 


diag(l,0+) 




o{2) 
0^ 


(1,1,0,0,0,0,0,0,0) 


{diag(l,-l,-l),diag(-l. 


1,-1)} 



(1,0,0,0,0,0,0,0,0) diag(l,0+) Udiag(-1,0_) 
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Table 4. Correspondence between Trebin and Longa's polar states and isotropy type strata for the 
group 0(3); the Graphical Representation (G.R.) of the typical point is the same as in |^: the straight 
lines denote the non degenerate eigenvector directions of Q, whilst the arrows represent the orientation 
of the polar vector P with respect to the eigenvectors of Q. The state Fij2 of [211 j for which the angle 
between the polar vector P and the unique non-degenerate eigenvector of Q is different from and 
7r/2, actually belongs to the stratum S^'^K Therefore, it is written in parenthesis, since it is equivalent 
(from the symmetry point of view) to the state -Fsg. 



Stratum 



isotropy subgroup 



G. R. of the typical point Phase in [2^ 



{e} 




ch 



{diag(±l,l,l)} 




{diag(l,±l,±l)} 



'B\\ 



s 



(2) 



{diag(l,0), Og 0(2)} 



S 



(2) 
B 



{diag(±l,±l,±l)} 



S« {diag(±l,0), O G 0(2)} 
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Figure 1. Possible phase transitions between bordering strata are connected by continuous sequences 
of one or more arrows. For each possible structural phase (isotropy type stratum), a "typical" point in 
the order parameter space and the corresponding isotropy subgroup of 0(3) are indicated. 
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